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Abstract: The antenna subtraction formalism allows to calculate QCD corrections to jet 
observables. Within this formalism, the subtraction terms are constructed using antenna 
functions describing all unresolved radiation between a pair of hard radiator partons. In 
this paper, we focus on the subtraction terms for double real radiation contributions to 
jet observables in hadron-hadron collisions evaluated at NNLO. An essential ingredient to 
these subtraction terms are the four-parton antenna functions with both radiators in the 
initial state. We outline the construction of the double real subtraction terms, classify all 
relevant antenna functions and describe their integration over the relevant antenna phase 
space. For the initial-initial antenna functions with two quark flavours, we derive the phase 
space master integrals and obtain the integrated antennae. 
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1. Introduction 

Jet production processes constitute an important tool for precision studies due to their large 
cross sections at high energy colliders. Reliable theoretical predictions for these observables 
require the calculation of at least the next-to-leading order QCD corrections. For these 
observables, the inclusive cross section with two incoming hadrons H\,H2 can be written 

as 

da = E/ d6d6/a/i(6)/ b /2(6)d(7 ab (6 H U (, 2 H 2 ) , (1.1) 

a,b 

where £i and £2 are the momentum fractions of the partons of species a and b in the incom- 
ing hadrons, / being the corresponding parton distribution functions and dcr^^i-Hi, £2-^2) 
is the parton-level scattering cross section for incoming partons a and b. 

The partonic cross section da a b has a perturbative expansion in the strong coupling a s 
such that theoretical predictions for a hadronic process at a given order in a s are obtained 
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when all partonic channels contributing to that order of the partonic cross section are 
summed and convoluted with the appropriate parton distribution functions as in eq. (1.1). 

In general, beyond the leading order, each partonic channel contains both ultraviolet 
and infrared (soft and collinear) divergences. The ultraviolet poles are removed by renor- 
malisation in each channel. Collinear poles originating from the radiation of initial state 
partons are cancelled by mass factorisation counterterms and absorbed in the parton dis- 
tribution functions. The remaining soft and collinear poles cancel among each other when 
all partonic contributions are summed over [1]. As jet observables depend in a nontriv- 
ial manner on the experimental criteria used to define them, they can only be calculated 
numerically. The computation of hadronic observables including higher order corrections 
therefore requires a systematic procedure to cancel infrared singularities among different 
partonic channels before any numerical computation of the observable can be performed. 

For the task of next-to-leading order (NLO) calculations, the infrared divergences 
present in real radiation contributions can be systematically extracted by process-indepen- 
dent procedures, called subtraction methods. The purpose of any subtraction method at 
NLO is to provide a subtraction term which has the same singular behaviour as the real 
radiation squared matrix element and is sufficiently simple to be integrated analytically 
over the radiation phase space which has been factorised from the (to + l)-particle phase 
space. The actual form of this subtraction term depends on the subtraction formalism used. 
Several successful subtraction formalisms have been proposed in the literature [2-6]. Most 
notably, the FKS [3] subtraction by Frixione, Kunszt and Signer and the dipole formalism 
of Catani and Seymour [2] have been implemented in an automated way, the former in [7], 
the latter in [8-13]. The major challenge for NLO calculations is the computation of one- 
loop amplitudes for multiparticle processes. The evaluations of 2 — > 4 processes at the 
next-to-leading order represent the current frontier [14-17]. 

Nevertheless, for some hadronic processes, in particular 2 — > 1 or 2 — > 2 scattering 
processes such as Drell-Yan, Higgs production, dijet production, vector-boson plus jet, 
vector-boson pair production or heavy quark pair production, the accuracy of the next-to- 
leading order predictions is not sufficient to match the anticipated experimental accuracy, 
expected to reach the order of a few percent or better. Accurate precision studies enabling 
the extraction of fundamental parameters of the theory will require that the theoretical 
predictions have the same precision. Those need therefore to be evaluated up to the next- 
to-next-to-leading order (NNLO) in perturbative QCD . 

The calculation of observables with m jets in addition to other objects (like for example 
vector bosons) at the NNLO requires three distinct contributions: the double real radiation 
d&NNLO w ith (to + 2) final state partons, the mixed real-virtual radiation da^ l NLO with 
(to + 1) final state partons and the 2-loop virtual radiation da^ NLO , with to final state 
partons. Those build the NNLO cross section which is given by 



The individual contributions in the to-, (to + 1)- and (to + 2)-parton final states are all 




(1.2) 
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separately infrared divergent. After renormalisation and factorisation, their sum is finite, 
though. For most massless jet observables of phenomenological interest, the two-loop 
matrix elements have been computed some time ago, while the one-loop matrix elements 
are usually known from calculations of NLO corrections to (m + l)-jet production [18] . 
The one-loop and two-loop matrix-elements contain explicit infrared divergences from the 
loop integration. Those cancel with divergences which are implicit in the real radiation for 
(m + 1)- and (to + 2)-parton processes. These real radiation divergences become explicit 
only once the phase space integration is carried out. The main issue of these calculations 
is therefore to find a method to extract and cancel the infrared divergences among these 
three contributions in order to finally evaluate numerically the finite remainders to obtain 
the NNLO contribution to the cross section. 

Like at NLO, a subtraction formalism is needed in order to extract the infrared diver- 
gences from these contributions. At parton level, the general form of the cross section for 
an m-particle final state at NNLO including subtraction terms is given by [19]: 



daNNLO = / { d °NNLO ~ ^NNLo) + / ^NNLO 

[ (a- V ^ -A-VS' 1 \ _l_ f r]~VS,l f ,.MF,1 

T / { aa NNLO aa NNLOj + I ® N N LO I ® NNLO 

+ / da NNLO+ / d& NNLO- ( L3 ) 

Here, da NNLO denotes the subtraction term for the (to + 2)-parton final state which 

behaves like the double real radiation contribution da^ NLO in all singular limits. Likewise, 
V S 1 

da N] j LO is the one-loop virtual subtraction term coinciding with the one-loop (to + 1)- 
particle contribution da^ LO in all singular limits. The two-loop correction to the to- 
parton final state is denoted by da^ NL0 . In addition, when there are partons in the 
initial state, there are two mass factorisation contributions, da^J^ and da^^ , for the 
(to + 1)- and TO-particle final states respectively. Like at the next-to- leading order level, 
the subtraction terms are needed in their unintegrated as well as in their integrated forms. 

There have been several approaches to build a general subtraction scheme at NNLO 
[19-29]. Another subtraction scheme, the (/^-subtraction formalism has been proposed 
in [30] . It has been applied to evaluate observables related to processes with colourless high 
mass final states [30-33]. In addition, there is a completely independent approach called 
sector decomposition which relies on a systematic expansion of the integrals in distributions 
followed by a purely numerical integration. It has been developed for virtual [34-36] and 
real radiation [37-40] corrections at NNLO, and applied to several observables already [41- 
44]. 

We will follow the NNLO antenna subtraction method which was derived in [19] for 
decays of a colourless initial state into massless final state partons. This formalism has 
been applied in the computation of NNLO corrections to three-jet production in electron- 
positron annihilation [45-48] and related event shapes [49-53], which were subsequently 
used in precision determinations of the strong coupling constant [54-58] . 
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For processes with initial-state partons, the antenna subtraction formalism has been 
so far fully worked out only to NLO in [59]. It has been extended to NNLO for processes 
involving one initial state parton relevant for electron-proton scattering in [60] while an 
extension of the formalism to include two initial state hadrons at NNLO is under con- 
struction [61,62]. An essential step towards this aim is performed in [62] where an explicit 
derivation of the subtraction terms needed for the double real contributions to the six- 
gluon process is presented. The general structure of the unintegrated subtraction terms 
relevant for the double real contributions to any hadronic observables evaluated at NNLO 
is presented there as well. 

In this paper, we will focus on the integrated form of the subtraction term relevant for 
double real radiation for processes involving two partons in the initial state. In the antenna 
subtraction formalism, the subtraction terms are built with so-called antenna functions. 
The latter describe all unresolved partonic radiation off a hard pair of colour-ordered 
partons, the radiators. Depending on where the two hard radiators are located, three 
cases need to be distinguished: both radiators are in the final state (final- final) , only one 
radiator parton is in the initial state (initial-final) or both radiator partons are in the initial 
state (initial- initial) . The subtraction terms and the antenna functions building them are 
separated corresponding to these three cases. In the most general hadronic process (two 
partons in the initial state, two or more partons in the final state), all three configurations 
have to be taken into account. 

As discussed in [19,60,62] the subtraction terms are separated according to the colour 
connection of the unresolved partons. In this paper, we will specialise on the subtraction 
term for the case where two unresolved (soft or collinear) partons are colour connected to 
the two incoming partons. This is indeed the only case where new ingredients, namely 
the four-parton initial-initial antennae are needed in unintegrated as well as in integrated 
forms. On a longer term we are aiming to evaluate the whole set of integrated four-parton 
initial-initial antennae, those are universal building blocks for the subtraction terms for 
any hadronic process evaluated at NNLO. 

In a first step towards this aim, in this paper, we have focused on the crossings of 
two partons from a set of three 4-parton final-final antennae involving two quark flavours. 
More precisely, the paper will be organised as follows: In Section 2, we present the gen- 
eral formulae for the subtraction terms related to double real radiation for initial-initial 
configurations while the colour-connected case is treated explicitly in Section 3. Section 4 
establishes a list of all non-identical initial-initial four-parton antenna functions relevant 
to construct the subtraction term for the double real radiation off two initial-state par- 
tons. In Section 5, the phase space mapping appropriate for initial-initial configurations is 
presented. 

Finally, Section 6 contains our results for the integrated initial-initial antenna functions 
with two quark flavours and Section 7 our conclusions. Since the results are lengthy, we 
show only the leading pole terms in the manuscript, and attach the complete results as a 
Mathematica file. 
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2. Antenna subtraction for double-real radiation at NNLO in the initial- 
initial configuration 



Antenna subtraction has been derived explicitly for final-final and initial-final configura- 
tions at NLO and NNLO in [19] and in [59, 60] respectively. For the initial-initial case, it 
has been derived at NLO in [59] and is under construction at NNLO. 

At NLO, the subtraction term is introduced to extract and cancel the infrared di- 
vergences present in the real contributions. The general forms of the subtraction terms 
required at NLO in any of the three configurations (final-final, initial-final and initial- 
initial) have been given in [19,59] and summarised in [62]. At this order, only tree-level 
three-parton antenna functions involving one unresolved parton are needed to build the 
subtraction terms. Those functions are usually denoted by X®j k , A? jfc , ■ in the three 
configurations. Their definitions will be recalled in Section 4. 

At NNLO, two contributions, double real and mixed real-virtual, require the intro- 
duction of a subtraction term. For final-final and initial-final configurations both types 
of subtraction terms have been constructed; at this order tree-level four-particle antennae 
involving two unresolved partons and one-loop three-parton antenna functions are needed 
respectively. Those functions have been derived and integrated over the corresponding 
factorised phase space in [19] and [60]. For initial- initial configurations, the general struc- 
ture of the unintegrated subtraction terms relevant for the double real contributions to 
pp — > m jets is presented in [62], we will recall it in this section. 

The double real radiation contribution to the m-jet cross section in pp collisions reads 

&°NNLO = N d ®m+2(kl, • • • , k m+2 ;p l ,p 2 )—^— 

m+2 

x \M m+2 (pi,p 2 ;k u ...,k m+2 )\ 2 4™ +2) (k 1 ,...,k m+2 ;pi,p 2 ) ■ (2.1) 

In this equation, \M m + 2 (pi,p 2 ; k\, . . . , k m+2 )\ 2 stands for the colour-ordered 2 — > m + 2 
matrix-element squared. The symmetry factor S m + 2 accounts for identical partons in the 
final state. The normalisation factor J\f includes all QCD-independent factors as well as 
the dependence on the renormalised QCD coupling constant a s . Ylm+2 denotes the sum 
over all configurations with m + 2 partons. The initial state momenta are labelled as usual 
as pi and p 2 whereas the m + 2 momenta in the final state are labeled k\ , . . . , k m+2 . d<3? m+2 
is the 2 — > m + 2 particle phase space 

d$ m+2 (£;i, . . ■ ,k m+2 ; Pl ,p 2 ) = [dfci] . . . [dk m+2 ] (2ii) d 5 d {p 1 + p 2 - h - ... - k m+2 ) (2.2) 

where we have introduced the abbreviation [dfci] = (d d_1 /ci)/(2E , i(27r) d_1 ). The jet func- 
tion jln l+2 \ki, . . . , k m+2 ;pi,p 2 ) ensures that out of (m + 2) final state partons, an observ- 
able with m jets is built. The incoming parton momenta pi,p 2 serve as reference directions 
to define transverse momenta and rapidities of the jets. 

The double real radiation contribution given in eq. (2.1) can become singular if ei- 
ther one or two final state partons are unresolved (soft or collinear). Consequently, when 
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constructing the corresponding subtraction term dafj NLO in eq. (1.3), which shall cor- 
rectly subtract all those single and double unresolved singularities we must distinguish the 
following configurations according to the colour connection of the unresolved partons: 

(a) One unresolved parton but the experimental observable selects only m jets. 

(b) Two colour-connected unresolved partons (colour-connected). 

(c) Two unresolved partons that are not colour connected but share a common radiator 
(almost colour- unconnected) . 

(d) Two unresolved partons that are well separated from each other in the colour chain 
(colour- unconnected) . 

(e) Large angle soft gluon radiation. 

This separation among subtraction contributions according to colour connection is 
valid in final-final, initial-final or initial-initial configurations and in any of those the sub- 
traction formulae have a characteristic structure in terms of required antenna functions. 
This antenna structure has been derived for the final-final and initial-final cases in [19,60] 
and in [62] for the initial-initial case. 

In here, we focus only on the initial-initial case with the kinematical situation where 
two unresolved partons are colour-connected to the two incoming partons. This is the only 
case where new ingredients are needed, namely the initial-initial four-parton antenna func- 
tions denoted by Xuj^. The four-particle initial-initial antenna functions will be defined 
explicitly in Section 4. 

3. Subtraction terms for two colour-connected unresolved partons in the 
initial-initial configuration 

When two unresolved partons j and k are adjacent and colour-connected to two initial- 
state partons, the subtraction term related to the double real contribution dajy NLO given 
in eq. (2.1) reads : 

d& NNLO = N d $m+2(h, • • • , fc m+2 1 i?l , f>2 ) 

^2 bm + 2 

X ^2^2 ( X il,jk ~ X l,jk X iL,K ~ X i,kj X Il,j) 
il jk 

x \M m (pf,pi; ki, . . . , fcj, fcj, ... , k m+2 )\ 2 J^Hh, ■ ■ ■ , h, h, ■ ■ ■ , k m +2;Pi,Pi) , (3.1) 

where the sum runs over all colour- adjacent pairs j, k and implies that the hard momenta 
i,l are chosen accordingly. X^ d - k denotes a four-particle tree-level initial-initial antenna 
function defined explicitly in Section 4. By construction those contain all colour connected 
double unresolved limits of the 2 — > m + 2 parton matrix element associated with partons 
j and k being unresolved between radiators i and I. 
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Figure 1: Illustration of NNLO antenna factorisation representing the factorisation of both the 
squared matrix elements and the (m + 2)-particle phase space when the unresolved particles j and 
k are colour connected between two initial state radiators i and I. 

However this antenna can also become singular in single unresolved limits associated 
with j or k only where it does not coincide with limits of the matrix element. To ensure 
that this subtraction term is only active in the double unresolved limits of the real matrix 
elements squared we remove these single unresolved limits of the four-particle antennae. 
Those limits are products of two tree-level three-particle antennae, namely products of an 
initial- final antenna X a ^ c and an initial-initial antenna X^c- I n these antennae we have 
replaced the original hard radiators with new particles, I and L. When both radiators 
are in the initial state as it is the case here, Pj = XiPi, p^ = x\p\. The product of these 
three-particle antenna functions in do^^o ^ nen subtracts the single unresolved limits of 
the associated four-particle antenna. 

The 2 — >■ m matrix element \M m {pf,Pi', k\, . . . , fcj, k\, . . . , k m+2 )\ 2 is evaluated with 
new on-shell momenta which are Lorentz boosted as a result of the mapping required to 
ensure factorisation of matrix element and phase space. The NNLO initial-initial mappings 
have been discussed in [59] and will be recalled in this paper in Section 5. 

Using a factorised form for both matrix element and phase space as shown in figure 
1, one is able to obtain an integrated form of the subtraction term with m partons in the 
final state which can be combined with the virtual two-loop contributions having also m 
final state partons as defined in eq. (1.3). 

More explicitly, the factorisation of the phase space with m + 2 final state particles 
denoted by k±, . . . , k m+ 2, reads, 

d$ m+2 (/ci, . . . , k m+2 ;pi,P2) = d$ m (fci, ...,ki,k h ..., k m+2 ; x 1 p 1 ,x 2 p 2 ) 

xS(xi — x\) 5(x 2 — x 2 ) [dkj] [dfefc] dx\dx 2 . (3.2) 

Using this factorised form of the phase space we can rewrite the integrated (colour- 
connected) subtraction term involving only the four-parton antennae in the form, 

\M m \ 2 Jt } d$ m J [dkj] [dk k ] 6{ Xl - Xl ) 5(x 2 - x 2 ) Xl jk d Xl dx 2 . (3.3) 

Moreover, the integrated antennae are defined as the antenna functions integrated over 
the antenna phase space as defined in eq. (3.2) including a normalisation factor to account 
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for powers of the QCD coupling constant by, 



1 



[dkj] [dk k ] Xi xi S(xi - Xi) 5(xi - xi) X iljk , (3.4) 



[C(6)] 2 



where C(e) is given by, 



C(e) = (4vr) 




(3.5) 



These integrations are performed analytically in d dimensions to make the infrared 
singularities explicit. The integrated initial-initial antennae are presently unknown and 
the calculation of a sub-set of those involving two quark flavours is the subject of Section 
6 and constitutes our main result in this paper. A list of all non-identical initial-initial 
four-particle antennae will be presented in Section 4. 

4. Initial-initial antenna functions 

In this section we shall recall how the antenna functions are defined in any of the configu- 
rations and how they enter in the construction of the subtraction terms. We also present 
a list of all non-identical four-parton initial-initial antenna functions. 

4.1 Definitions 

In the antenna formalism, in any of the three configurations (final-final, initial-final or 
initial-initial) the subtraction terms are constructed from products of antenna functions 
with reduced matrix elements (with fewer final state partons than the original matrix 
element). The integrated subtraction terms are obtained after an integration over a phase 
space which is factorised into an antenna phase space (involving all unresolved partons 
and the two radiators of the antenna) multiplied with a reduced phase space (where the 
momenta of radiators and unresolved radiation are replaced by two redefined momenta). 
These redefined momenta can be in the initial state, if the corresponding radiator momenta 
were in the initial state as we saw in Section 3 in eq. (3.1) for the subtraction term v^nlo- 
An antenna function is determined by the external states it contains and the pair of 
hard partons it collapses to in the unresolved limits. In general we denote the antenna 
function as X. For antennae that collapse onto a hard quark- antiquark pair, X = A for 
qgq. Similarly, for a quark-gluon antenna, we have X = D for qgg and X = E for qq'q' 
final states. Finally, we characterise the gluon-gluon antennae as X = F for ggg, X = G 
for gqq final states. 

At NLO, we only need to consider tree level three-particle antennae involving only one 
unresolved parton. At NNLO we will need four-particle antennae involving two unresolved 
partons and one-loop three-particle antennae. 

In all cases the antenna functions are derived from physical matrix elements associated 
to the decay of a colourless particle into partons: the quark-antiquark antenna functions 
are derived from 7* — > qq + (partons) [63], the quark-gluon antenna functions from x ~ > 
g + (partons) [64] and the gluon-gluon antenna functions from H — > (partons) [65]. The 
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tree-level antenna functions are obtained by normalising the three- and four-parton tree- 
level colour sub-amplitudes squared to that of the basic two-parton process: The final-final 
three- and four-particle antennae are respectively defined by: 



Sj.-j 



ijk,IK 



X, 



\M% k \* 
\M°j I 2 

\Ml 



l IK\ 



ijkl 



Hjkl,IL 



ijkl I 



\M 



(4.1) 



IL\ 



where S denotes the symmetry factor associated with the antenna, which accounts both 
for potential identical particle symmetries and for the presence of more than one antenna 
in the basic two-parton process. It is chosen such that the antenna function reproduces 
the unresolved limits of a matrix element with identified particles. 

The initial-final tree-level three- and four-parton antennae denoted by X^ - k and X^ - kl 
are in principle obtained by crossing one-parton to the initial state starting from the cor- 
responding final-final antennae. However this crossing might be ambiguous as was first 
noticed in [59] for the quark-gluon type antenna. Depending on the unresolved limit con- 
sidered, the pair of hard partons it collapses to may be different. 

The initial-initial tree-level three- and four-parton antennae denoted by X® k t and X^ - k 
are obtained by crossing two partons to the initial state, starting from the final-final an- 
tennae. This crossing procedure, unlike in the initial-final case, is free of ambiguity as the 
pair of hard partons the initial-initial antenna collapses to is always uniquely defined. 

More explicitly, the three-parton initial-initial antenna function Xf k ■ is defined as 

M (Xj,k] ^ 

Q 

°ijk,IK~ 



Y° 

A ik,j 



M\I,K 



(4.2) 



where i denotes that particle i is crossed to the initial state. Therefore we can see that the 
initial-initial antenna is connected to the final-final antenna (where all coloured particles 
are outgoing) by 

X? kJ = (-l) A FXl(- Pi , Pj ,-p k ) (4.3) 

with Ap the difference in the number of fermion crossings between the three-particle and 
the two-particle subamplitude. 

The tree-level four-particle initial-initial antenna X^ - k is defined as 

M(i,j,k,l) 



X 



il,jk 



s, 



ijkl, I L' 



M(I,L) 



and the relation to the final-final antenna is 



A,- 



ljk = {-ir F x u i {-p i ,p j ,Pk-pi) 



(4.4) 



(4.5) 



The four-particle tree-level antenna functions are not determined by the species of the 
particles alone but also by the colour-connection. We distinguish leading-colour antennae, 
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denoted by letters without tilde, where the particles are colour-connected in the order they 
are listed and subleading colour antennae, denoted by letters with tilde, where the gluons 
are photon-like. This notation has been used in [19, 60] and we will further use it in this 
section in the tables below. 

The unresolved limits of the initial-initial antennae can be obtained from those of the 
final-final antennae by crossing. The crossing of the triple-collinear splitting functions is 
explained in [66]. Those limits will be reported elsewhere. 

4.2 Catalogue of tree-level four-particle initial-initial antenna functions 

Any two particles of a four-particle final-final antenna can be crossed to the initial state 
to obtain an initial-initial antenna; therefore one final-final four-particle antenna gives rise 
to six initial-initial antennae. Due to symmetries, at most four of these initial-initial an- 
tennae are different. The independent crossings are listed below. To make the colour 
connection clear, in this list we write out the arguments of the antennae explicitly, i.e we 
write X® (i,j,k,lj (where i denotes an incoming particle) instead of X® k ^. Some initial- 
initial antennae are free of singular limits. These finite antennae are not needed for the 
construction of subtraction terms, but their integrated form could be needed for cross- 
checks. They are marked with ** below. 



quark-antiquark antennae 

A° Al (q, g, g, q) , A\ (q, g, g, q) , A° 4 (q, g, g, q) , A\ (q, g, g, q) 

A° 4 A\ (q, g, g, q) , A% (q, g, g, f ) , A° 4 (q, g, g, q) 

B° 4 B 4 {q,q',q',q), B° 4 (q, q> ,q> ,?) , B° 4 (q,q> ,§> ,q)** 

c° c l (?> q,q,q),c2 (?, q, q, q) ,c%(q, % q,q)**,c2(q, q, q, 



A\{q,g,g,q) is symmetric under the interchange of the two photon-like gluons. The 
nonidentical-flavour antenna B® is separately symmetric under interchange of q' with q' 
and of q with q. 



quark-gluon antennae 





D° 4 (q, 


9,9,9), Dl(q,g,g,g), D° A (q, 


9,9,9), E>l(q,g,g,g) 




El 


E° 4 (q, 


q',q',g),El(q,q',q',g),El( 


q,q',q',g), El(q,q',q', 


9) 


El 


El(q, 


q',q',g),El(q,q',q',g),El( 


q,q',q',g), El{q,q',q', 


9) 



Due to the cyclic colour connection, D A is symmetric under interchange of the second and 
fourth gluon. 
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gluon-gluon antennae 

F 4 ° F°(g,g,g,g), F°(g,g,g,g) 

Gl G\ (g, q, q, g) , G% (g, q, q, g) , G% (g, q, q, g) , G\ (g, q, q, g) 
G\ G% (g, q, q, g) , G\ (g, q, q, g) , & A (g, q, q, g) 
Hi H2(q,%q',q>), H° A (q,q,q',q>) 



F® is symmetric under cyclic interchange of its arguments. G4 is symmetric under the 
interchange of the two gluons as well as under the interchange of q with q. has three 
symmetries, q f-> q, q' -H- q' and the flavour renaming q q' , q -H- q' . 

5. Phase space factorisation and mappings 

The construction of subtraction terms requires a mapping from the original set of momenta 
onto a reduced set. The mapping interpolates between the different soft and collinear lim- 
its which the subtraction term regulates. An appropriate mapping for the initial-initial 
case, both for single and double unresolved configurations, has been discussed in [59]. By 
requiring momentum conservation and phase space factorisation, the phase space mapping 
is strongly constrained. The remapping of initial state momenta can only be a rescaling, 
since any transversal component would spoil the phase space factorisation. For two unre- 
solved partons j and k, a complete factorisation of the phase space into a convolution of 
an m-particle phase space depending on redefined momenta only and the phase space of 
the unresolved partons j and k can be achieved with a Lorentz boost. This boost maps 
the momentum q = p\ + P2 — kj — k^ , with q 2 > and pi,P2 being the momenta of the 
hard emitters, into the momentum q = x±p± + X2P2 , where x\ and X2 are fixed in terms 
of the invariants as follows: 

( S\2 - Sj2 - Sk2 S\2 - S\j - S\k ~ Sj 2 ~ Sfc2 + Sjk 

V s 12 S12 - Sij - Sifc 

_ ( S\2 — Slj — Sjk Su - Sij — Sik — Sj2 — Sk2 + Sjk 

V s 12 -512 - S j2 ~ S k2 

These two definitions guarantee the overall momentum conservation in the mapped mo- 
menta and the correct soft and collinear behaviours . The two momentum fractions x\ and 
X2 satisfy the following limits in double unresolved configurations: 

1. j and k soft: x\ — >■ 1, X2 — > 1, 

2. j soft and kj~ = z\py. x\ — >■ 1 — z±, X2 — > 1, 

3. kj = z\p\ and k k = z 2 p2- X\ ->■ 1 — z\ , x 2 ->■ 1 - z 2 , 

4. kj = zipi, k k = Z2P\: X\ ->■ 1 — z\ — z 2 , x 2 ->■ 1, 




and all the limits obtained from the ones above by the exchange of p\ with p2 and of kj 
with k k . The construction of NNLO antenna subtraction terms also requires that all single 
unresolved limits of the four-parton initial-initial antenna functions Xuj k , with radiators 
i and I, have to be subtracted, such that the resulting subtraction term is active only in 
its double unresolved limits. A systematic subtraction of these single unresolved limits by 
products of two three-parton antenna functions can be performed only if the NNLO phase 
space mapping turns into an NLO phase space mapping in its single unresolved limits. A 
detailed discussion of the corresponding translation between these two momentum map- 
pings can be found in [59]. 

The factorisation of the (m + 2)-parton phase space into an m-parton phase space and 
an antenna phase space involving the unresolved partons j and k given in eq. (3.2) can 
equivalently be written as 

d$ m+2 0i, . . • , k m+2 ;pi,P2) = d$ m (fci, ...,k u ki,..., k m+2 ; x 1 p 1 ,x 2 p2) 

x J 5{q 2 - x lX2 s 12 )6(2(x 2 p 2 ~ xipi).q) (5.2) 
x [dkj] [dkk] dxi dx 2 , 

where J is the Jacobian factor defined by 

J = Sl2 (Xl(si 2 - Slj - si k ) + x 2 (si 2 - s 2j - s 2k )) . 

This phase space parametrization can also be used to give an equivalent definition of 
the integrated initial-initial antennae first given in eq. (3.4). Those are given as, 

Xl lk ( Xl ,x 2 ,e) = J [dkj][dk k ] J Xl x 2 SiCJSWX^, (5.3) 

where 

Ci = q 2 - x 1 x 2 S12 , 

C 2 = 2(x 2 p 2 - xipi).q . (5.4) 

We shall use this definition of the integrated antennae j k {xi, x 2 , e) to compute them 
in Section 6. 



6. Integration of the four-parton initial-initial antennae 

In the first part of this section we describe how the initial-initial four-parton antennae are 
integrated over the antenna phase space while in the second part of this section we restrict 
ourselves to the evaluation of integrated antennae involving two quark flavours. 

6.1 Calculational method 

The initial-initial antenna functions have the scattering kinematics 

Pi + P2 -» kj + k k + q , 
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where q is the momentum of the outgoing colourless particle. The momenta satisfy: 



Pi = p\ = 0, fcf = A| = 0, q 2 = q 2 = x\ x 2 S12 • 



The four-parton initial-initial antennae defined in Section 4 need to be integrated over 
the phase space of the unresolved partons j and k. This integration yields a result which 
depends only on s\ 2 , x\ and x 2 . From dimensional counting, one can immediately conclude 
that the dependence on s\ 2 is only multiplicative, according to the mass dimension of the 
integral. 

The initial-initial antenna phase space integrals are derived from squared matrix elements 
and can be represented by forward scattering diagrams as in the following figure: 



The two delta functions in eq. (5.2) can be represented as mass-shell conditions of fake 
particles and are shown in the previous picture as a thick solid line (representing a massive 
particle with mass M 2 = x\ x 2 s\ 2 ) and a dashed line (representing the other constraint) . 
This allows us to use the optical theorem to transform the initial-initial antenna phase space 
integrals into cut two-loop integrals and therefore use the methods developed for multi- 
loop calculations [67-69] . Up to eight-propagator integrals with four cut propagators are 
generated in this way . Using the reduction techniques, the calculation of the integrated 
antennae can be related to the evaluation of a reduced set of master integrals . For the 
complete set of non-identical initial-initial four-parton antennae tabulated in Section 4, we 
find 32 such integrals, obtained using integration-by-part (IBP, [70]) and Lorentz invariance 
(LI, [71]) identities, following the Laporta algorithm [72]. A private implementation as well 
as a public one [73] have been used. 

For those four-parton antenna functions there are 13 different propagators, including the 
four that are cut in the phase space integration {Dj, Dk, £^fci2, Djki2s)'- 




Dji = {pi - kj) 2 , 
D k i = {pi - k k ) 2 , 
Dj2 = (P2 ~ kj) 2 , 
D k 2 = {P2 ~ k k ) 2 , 
Djk = {kj + k k ) 2 , 
Djki = {pi - kj - k k ) 2 
D jk 2 = {P2 - kj -k k ) 2 
Dji2 = {pi +P2- kj) 2 
D k 12 = {Pi +P2- kk) 2 
Dj = k 2 , 
D k = k 2 , 
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Djkl2 = (Pi +P2- kj - k k ) - X1X2S12 
? ifcl23 = (P3 + Pi + P2 ~ kj 



Djkl23 = (P3+P1 +P2- kj - k k ) 2 , (6.1) 



where p% = X2P2 — xiPi- To perform the reduction to master integrals, we drop any 
integral where Dj, D k , Dj k i 2 , Djki23 are not in the denominator and impose momentum 
conservation. The integrands of the 32 master integrals found can all be written as rational 
polynomials of the denominators above. 

These master integrals are calculated using either the method of differential equations 
or by a direct evaluation of the phase space integrals in terms of hyper geometric functions. 
The simplest master integral is a two loop bubble with all the internal lines cut, it is ob- 
tained from eq. (5.3) by replacing the jacobian J and the antenna with unity: 



= h{xi,x 2 ) = J d d q d d kj d d k k 5 d (pi + p 2 - q - kj - k k ) x 

5 + (kj) 5 + (k 2 k ) 6+ (q 2 - M 2 ) 6(2 (x 2P 2 - x lPl ).q) . (6.2) 
This integral can actually be expressed as a hyper geometric function: 



. (4tt)- 4+2£ r(l - 
* 12 r(2 - 2e) r(2 - 2e) 



2e \ I *■ \ x c ; -en „ \l-2e/i , „ \-e 



x -z(i- Xl y-^(i + Xl ) 



£ (1 - x 2 ) 1 - 2£ (l + x 2 )- e {x l + x 2 )- 1+2£ (6.3) 

(i- xi )(i-x 2 y 



2 F 1 e,l-e,2-2e, 



(l + si)(l + x 2 ) 



it has been checked against the master integral I[0] appearing in the calculation of the 
gauge-boson rapidity distribution at NNLO in [68], the notation translates as u = x 2 /x\, 

Z = X\X2- 

The set of master integrals which we denote by Ii(x%, x 2 , e) are functions of xi, X2 and e. 
We begin by factoring out the leading behavior of the master integrals h(x\, x 2 , e) in the 
limits x\ — > 1 and x 2 — > 1, keeping the exact e-dependence: 

I t (x 1 ,X2,s) = (l-x 1 ) m i- 2e (l-X2) m2 - 2e T l (x 1 ,x 2 ,e). (6.4) 

The integers m\,m 2 are characteristic to each master integral. The functions J r i(xi,x 2 ,£) 
are regular at x± = 1, at x 2 = 1, and at x\ = x 2 =1 and can be calculated as Laurent 
series with, at most, second order poles in e. 

The integrated antennae given by X(xx,x 2 ,e) are linear combinations of these master 
integrals Ii(xi,x 2 , e), with coefficients containing poles in e, as well as in (1 — x\) and 
(1 — x 2 ). After the masters have been inserted into the integrated antennae, those take the 
form 

X(x!,x 2 , e) = (1 - x^- 1 - 2 ^! - x 2 )~ l - 2£ K(xi,x 2 ,e). (6.5) 



-14- 



where TZ(xi, X2, s) is regular at the boundaries x\ = 1, X2 = 1, and at x\ = X2 = 1. The 
e-expansion of the singular factors (1 — Xj)~ 1 ~ 2e is done in the form of distributions: 

(1 - x t y^ = -1 6(1 - Xi ) + V , ( 6 .6) 

is n! 

n 

with 

To evaluate the integrated antennae, we decompose the phase space into four regions 
depending on the values of x\ and xi- Those regions are given by: 

• X\ 7^ 1, X2 7^ 1, which we refer to as the hard region 

• a?i = 1, £2 7^ 1, and a?i 7^ 1, X2 = 1, referred to as collinear regions 

• X\ = 1, #2 = 1, which we denote the soft region. 

In the hard region (x\ 7^ 1, X2 7^ 1), harmonic polylogarithms of weight two appear in the 
O(e ) term of 7?-. Therefore, the e-expansion of the master integrals in the hard region is 
needed at least up to the order at which terms of transcendentality 1 two appear for the 
first time generally. 

In the collinear regions (x\ = 1 or X2 = 1), since the expansion in distributions (6.6) 
generates additional 1/e factors, the function 1Z is required up to 0(e) where harmonic 
polylogarithms of weight 3 appear. The masters evaluated in the collinear region need 
therefore to be expanded at least up to the order at which terms of transcendentality 3 
appear generally. 

Finally, in the soft region (x\ = X2 = 1), since the expansion of the distributions (6.6) 
generates additional 1/e 2 coefficients, the function 1Z is required up to 0(e 2 ) where poly- 
logarithms of weight 4 appear. The masters evaluated in the soft region need therefore to 
be expanded up to transcendentality 4 at least. 

6.2 Integrated antennae with two quark flavours 

In a first step towards the calculation of all integrated four-parton initial-initial antenna 
functions for the double real radiation case, in this paper, we focus on all the crossings of two 
partons from the following four-parton final-final antennae: B®(q,q' \q' ,q), E®(q,q',q',g) 
and H®(q, q, q' , q') defined in [19] . 

We found that the reduction involving only those initial-initial antennae leads to 12 
master integrals. Those without numerators are shown in Fig. 2 . We have performed the 
calculation of these integrated antennae with two choices of master integral bases differing 
by four master integrals. In the first basis, the definitions of the master integrals involved 
in the calculation are as follows: 

h = J [dkj] [dk k ] 5(C 1 )5(C 2 )(h ■ Pl ) , (6.8) 



1 We define the transcendentality of a product as the sum of the transcendentalities of its factors, GHPL's 
and HPL's have transcendentality equal to their weight, log 2 has transcendentality one, £(n) has transcen- 
dentality n. 
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h = J [dkj] [dk k ] 8{C{)8{C 2 ) , 

J L>kl2 

h = j [dkj] [dfc fc ]«5(C 1 )«5(C 2 )^-^, 
h = J [dk^ [dfc fc ] (5(Ci)(5(C 2 ) ~ ( ^' P3) . 



fcl2 



J 6 = / [d^tdfcfel^d)^)-^-, 

I 7 = [ [dkj] [dfc jfc ]«5(C 1 )«5(C 2 )-^-, 
J 8 = y [dfcj] [dfc fe ]<5(C\)<5(C 2 )-^-, 
/ 9 = / [dfc^[dfc fc ]$(Ci)5(C 2 )-^-, 

J u k\ 

I l0 = [ [dkj] [dfc fc ] <5(Ci)<5(C 2 ) 



7i4= / [d^Hd^]^)^)— - 4 

J ^12 J 

iis = y [dfcj] [dfc fc ]5(Ci)*(C 2 )-^ 



l jlDk2 

In the other choice of basis, the masters with scalar products in the numerator (Ji, ^3, J4, is) 
are replaced by alternative master integrals: 



f( = y [dfcj] [dfc fc ] 8{C{)8{C 2 )D 



i' 3 = / [dfcjHdfcfcWCiWCb) 



fcl2 

,jDfc12 (6.9) 

i' 4 = J [dkjUdkkWCMc*)^- 
i' 5 = j [dfc j ][dfc fc ]«y(c 1 )«y(c 2 )— -1— ■ 

J L) ik L) k\2 

The relation to the first basis is: 

7i = si 4 /2 - 2(,r+* 2 ^ 

2xf + 9x1^2 +xk- \2x\x\ + 2e(-2x? - 7xix 2 - x^ + lOxf s|) r ,„ 1 n . 

73 = 512X2 ^K^^ /2 (fU0) 



-512^1 



-xf + xix 2 — ~~ 4xfx^ + 8x1X3 
8(-l + 2e)(xi -x 2 )(xi +x 2 ) 

x\ — XlX 2 + 2x| + 4x^X2 — 6xix| 



+ 2e 



8(-l + 2e)(xi -x 2 )(xi +x 2 ) 



+ 512^(^1 + X 2 + 4xiX2)/ 4 
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Figure 2: Master integrals for the phase space integration of the tree-level initial-initial B4, Hi and 
£4 type antennae at NNLO. Thick solid and dashed lines refer to the conditions on the phase space 
integral implemented as auxiliary propagators. All the internal lines are massless except for the thick 
solid line. Only the integrals without numerators are shown in this picture. 



s 3 f 
s l2 e 



(-1 + Xl)X\(l + Xi)(-1 + X2)X%{1 + £ 2 )(-l + XiX 2 ) 



+ S 12 Xl 



-(xi 



2(-l + 2e)(xi -x 2 )(x 1 +x 2 ) 
x 2 ){xi + x 2 )(l + 2xl) 



-1 + 2e)(a?i - x 2 ){xi + x 2 ) 



+2e- 



Q T 2 j_ 922 9 4 1 924 

-1 + 2e)(xi - X2){xi + ac 2 ) 



X X , X!+X 2 , Xi 

h = -y J 3 + — j — 4 + 12 Y 6 



4 



Before we proceed with the details of calculating the masters, we present in Table 1 a sum- 
mary of which regions contribute to the crossings of the antenna functions B®(q, q' , q 1 ,q), 
E®(q,q',q',g) and H®{q,q,q' \q'). We then present results of the required masters in each 
of these regions. As explained in Section 4.2, not all of the six crossings of a final-final 
antenna are different. Labelling the final-final antenna functions as B®(lq, 3q' , 4(f, 2q), 
E2(lq,2q',3q',4g) and H%(lq,2q,3tf,4(?), the identical crossings are the following: 



H12 


= H34 






His 


= H14 


— H23 


— H24 


B13 


= Bu 


= B23 


= B24 


E\2 


= E13 






E24 


= E34 
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where 



Hl2 Hqq qtqt , 


(6.12) 


^13 = Hqq',qq' ' 


(6.13) 


-^i2 qq,q q i 


(6.14) 


^34 = B^ q/ qq, 


(6.15) 


Bl3 = B qq',q'q, 


(6.16) 


E\2 = E °qq',q'g , 


(6.17) 


Eu = Eq g ,q'q > 


(6.18) 


-^23 = Eq, ql q , g , 


(6.19) 




(6.20) 



As mentioned in Section 4.2, B34 is free of singular limits. It will therefore not be needed 
for the construction of subtraction terms and its integrated form is free of poles in e. It 
might be needed for checks of the integrated antennae, however, which is why it is included 
here. 



Antenna 


soft 


collinear x\ = l 


collinear X2=l 


hard 


H12 


no 


no 


no 


yes 


n 13 


no 


no 


no 


yes 


£12 


no 


yes 


no 


yes 


£u 


no 


yes 


no 


yes 


£23 


no 


no 


no 


yes 


£24 


no 


no 


no 


yes 


612 


yes 


yes 


yes 


yes 




no 


yes 


no 


yes 


£34 


no 


no 


no 


yes 



Table 1: Summary of the regions contributing to each of the independent crossings of the three 
antennae: B\, £4, H®. 



6.2.1 Master integrals 

In the following, we present the results for the masters in the hard, collinear and soft 
regions while restricting ourselves to those which are explicitly involved in the calculation 
of the integrated antennae with two quark flavours. 
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a) The hard region 

The master integrals denned above were computed in the hard region mainly with the 
differential equations technique [71,74,75]. The only masters that were calculated directly 
are I±, I 2 , and Since the dependence of the integrands on x\ and X2 is via the 

constraints, C\ and C 2 , as shown in eq. (5.4), we derive the differential equations for each 
master integral by employing the following operators at the integrand level: 

_^ = dC 1 J_ dC2^_ 

dx x dx 1 dCx dx x dC 2 ' y ' ' 

d dC x d dC 2 d . 

= 1 . (6.22) 

dx 2 dx 2 dC\ dx 2 dC 2 

The boundary conditions required for the solution of the differential equations are either 
obtained from self-consistency conditions on the integrals, or by explicit evaluation in the 
collinear or soft limits. The solution of the system of differential equations yields two- 
dimensional generalized harmonic polylogarithms (GHPL, [76,77]) of up to weight two, or 
products of weight one harmonic polylogarithms (HPL, [78]) of argument x\ or x 2 . The 
definition of the HPL and GHPL functions involved in the solution of the master integrals 
in the hard region is recalled below: 

H(l,a;) = -ln(l-x) 
H(0,x) = In x 
H(-l,x) = ln(l +x) 

The harmonic polylogarithms of higher weight are defined recursively (we group weights 
into vectors, b\, . . . ,b w = b): 

B(6 w ,x) = ^-\n w x (6.23) 
w\ 



while, if a = (a, b) / V 
with weight functions 



R(a,b,x)= [ dzf(a,z)H(b,z) 
Jo 



f(0,z) = - z , (6.24) 

This results in the derivative formula 

d 

—H(ai,b,x) = f(a!,x)H(b,x) 

The two-dimensional generalized harmonic polylogarithms are defined in a very similar 
fashion: 

G(0,y) = \ny , 
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G(l,y) = ln(l-y) , 
G(-l,y) = ln(l + y) , 
G(-z,y) = ln(l + ^ 



For weight w>lwe have 



G(O w ,y) = —\n w y , 
w\ 

G(a,y) = f V dy'g(a,y') G(b,y') 
Jo 



where 



y) 



1 



y-a 



and we have for the derivatives 



d_ 

dy 



G(a,y) = g(a,y)G(b,y) . 



(6.25) 

(6.26) 
(6.27) 

(6.28) 
(6.29) 



As can be seen above, two-dimensional generalized harmonic polylogarithms reduce to 
harmonic polylogarithms if their weights are only 0, 1, —1, with 



G(a, x) 



— H(a, x) if a contains an odd number of 1 
H(a, x) else. 



(6.30) 



In the results for the hard region below, the GHPL's will have x\ as their argument and 
0, ±1 and — X2 as their weights, the conversion to HPL's has been inserted where possible. 
An important check on our results in this region is to compare the soft and collinear limits 
of the hard result against the ones derived in the soft and collinear regions calculated by 
a direct integration of the phase space. We found agreement on all the powers of the e- 
expansion in the hard region using these two different methods. The results of the masters 
expanded to the needed order in e are given below. The common prefactor N? is defined 
as: 

(4tt)- 4+2£ 



r(2-2e) 



.31 




xk k -pi =I 1 = (s 12 ) 1 ' 



_ 2£ N r (l- Xl )- 2£ (l-x 2 
8(x 1 +x 2 ) 2 
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(xi - l)(x 2 - l) 2 (x 2 xi + 2x 1 + x 2 ) 

+ e |(xi - l)(x 2 xx + 2xi +x 2 )H(0,xi)(a; 2 - l) 2 
(5x1^2 + 3x2 + 3xiX2 — 7x 2 — 12xi) (x 2 — 1) 



1 x 



2 (x 2 X2 



2x\x 2 — 3x 2 x 2 — x 2 — 2xi) G(— X2,xi) 



+ (xi + l)(x 2 + l) 2 (x 2 xi - 2xi - x 2 )(H(-l, x\) + H(-l, x 2 ) 
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-H(0,x 2 )) + 4(x 1 +x 2 ) 2 log(2)} 

+ e 2 - l)(x 2 xi + 2xi + x 2 )(H(0, -1, x 2 ) - H(0, 0, xi) - H(0, 0, x 2 )) x 

(x 2 - l) 2 — -(xi — 1) (3xix 2 + x 2 + xix 2 - 5x 2 - 8xi) (x 2 - 1) 

+ - (xi — 1) (5xix 2 + 3x 2 + 3xix 2 — 7x 2 — 12xi) H(0, xi)(x 2 — 1) 

+ (-5x 2 x^ + 3x 2 + 8xxx 2 + 15x 2 x 2 + 7x 2 + 12xi) G(-x 2 , x\) 
+ (x 2 x 2 — x 2 — 2xix 2 — 3xfx2 — x 2 — 2xi) (2G(0, — x 2 , xi) 
+ 2G(-x 2 , -1, xi) + 2G(-x 2 , 0, xi) - 4G(-x 2 , -x 2 , xi) 
+ G(-x 2 , xi)(2H(-l, x 2 ) - 2H(0, x 2 ))) + (xi + l)(x 2 + 1) x 

(5xix 2 — 3x 2 — 3xix 2 — 7x 2 — 12xi) f—H(—l,xi) + -H(— l,x 2 ) 

-^H(0, x 2 )^ + {-x\x\ + x 2 + 2xix 2 + 6x 2 + 3x 2 x 2 + 12x x x 2 

+x 2 + 6x 2 + 2xi) (G(l,-x 2 ,xi) + H(-l,x 2 )H(l,xi) + H(l,-l,xi) 
+ H(l, — 1, x 2 )) + [x\x\ — x 2 — 2x\x\ — 6x 2 — 3x 2 x 2 — 12xix 2 
-x 2 - 6x1 ~ 2xi) (H(0, x 2 )H(l, xi) + log(2)H(l, xi) + H(l, 0, x 2 ) 
+ H(l,x 2 )log(2)) + (x! + x 2 ) 2 (-41og(2)H(0,x 1 ) -41og 2 (2) 
-4H(0, x 2 ) log(2)) + (xi + l)(x 2 + l) 2 (x 2 xi - 2xi - x 2 ) x 
(3G(-l,-x 2 ,xi) + H(0,xi)H(0,x 2 )-2H(-l,-l,xi) 
-2H(-l,-l,x 2 )-H(-l,0,xi) + 2H(-l,0,x 2 )-H(0,-l,xi) 
+H(-l,x 2 )(log(2)-H(0,xi)) + H(-l,xi)(-2H(-l,x 2 ) 

7T 2 

+2H(0,x 2 ) + log(2)) + y 

- 2(xi + x 2 ) (xix^ - 5x 2 - 6xi) log(2)} + 0(e 3 ) ] . (6.32) 



2e N r {l-x 1 )-*(l-x 2 )- 2e 
{12) (l-3e) (xi + x 2 ) 



+(xi - 1) (x 2 - 1) + e 1 2 log(2) (xi + x 2 ) + 2(xix 2 + 1) G(-x 2 , x x ) 

- (xi - 1) (x 2 - 1) H(0, xi) + (x! + 1) (x 2 + 1) (- H(-l, xi) 

- H(-l,x 2 ) + H (0,x 2 ))} 

+ e 2 {(xix 2 + l)(-2G(0,-x 2 ,xi) - 2G(-x 2 , -1, xi) - 2G(-x 2 , 0, x x ) 

+ 4G(-x 2 ,-x 2 ,X!) +G(-x 2 ,x 1 )(2H(0,x 2 ) -2H(-l,x 2 ))) 
+ (xi - 1) (x 2 - 1) (- H(0, -1, x 2 ) + H(0, 0, xi) + H(0, 0, x 2 ) ) 
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+ (x 2 xi+3xi + 3x 2 + l)(G(l,-x 2 ,xi) + H(-l,x 2 ) H(l,xi) 
+ H(l,-l,xi) + H(l,-l,x 2 )) + (xi + 1) (x 2 + 1) x 

(-3G(-l,-x 2 ,xi) - H(0,xi) H(0,x 2 ) + 2H(-1, -l,xi) + 2 H(-l, -1, x 2 ) 
+ H(-l,0,xi) -2H(-l,0,x 2 ) + H(0,-l,xi) + H(-l, x 2 ) (2 H(-l, xi) 



+ H(0,xi) -log(2)) + H(-l,xi)(-2H(0,x 2 ) - log(2) ) 



3~ 



+ (-x 2 x 1 -3xi-3x 2 -l)(H(0,x 2 ) H(l,xi) +log(2) H(l,xi) 
+ H(l,0,x 2 ) + H(l,x 2 ) log(2)) +(xi + x 2 ) (-21og(2) H(0,xi) 



-21og 2 (2) -2H(O,x 2 )log(2))}+0( £ 



(6.33) 



x(-(k rPl )(k rP3 )) =/ 3 = ( Sl2 ) i 8(xi + x 2 ) 2 X 

(xi - l)(5xi - 1)H(0, x 2 )x\ - 2xi (xi - x 2 )H(0, xi)H(0, x 2 )x\ 

- (xi - l)(5xi - x 2 )(x 2 - l)x 2 + xi(x 2 - 1) (-4x2 + x i x 2 + 2x 2 + %i) H(0, Xl) 
+ (-x? + 4x 2 xi - 6xi + l) H(0,xi)H(0,x 2 )x^ + (x x + l)(5xi + 1) x 

(-H(-l,xi)H(0,x 2 ) - H(-l,0,xi) - H(0,-l,x 2 ))x^ + (x x - l)(5xi - 1) x 
(2H(0, 1, x 2 ) - 2H(0, x 2 )H(l, x 2 ))x 2 2 + xi(12H(0, x 2 ) log(2) - 12H(0, x x ) log(2))x| 
+ xi(xi -x 2 )(-2G(0, -x 2 ,0,xi) -2G(0,-x 2 ,xi)H(0,x 2 ) 
-4G(-x 2 , 0, xi)H(0, x 2 ) + H(0, x 2 ) (2H(0, -1, xi) + 6H(0, 0, xi) - 4H(0, 1, xi) 

+— J + H(0, xi) (H(0, x 2 )(4H(l, xi) + 4H(1, x 2 )) + 2H(0, -1, x 2 ) 

-4H(0, l,x 2 ) + — J + 2H(0, -1, 0, xi) + 2H(0, -1,0, x 2 ) - 2H(0, 0, 0, x 2 ) 

+3C(3)) xl - (xi - l)(23xi - 5x 2 )(x 2 - l)x 2 - 2(5xi - x 2 )(xix 2 + 1) X 
G(-x 2 , xi)x 2 + (xi + l)(5xi - x 2 )(x 2 + 1)(H(-1, xi) + H(— 1, x 2 ))x 2 

+ (4x 2 X 2 — 6x 2 + x\x\ — 15X2X1 — 4xi + x \ + 3X2) H(0, X2)X2 

- 2(5xi - x 2 )(x\ + x 2 ) log(2)x 2 + ^7r 2 (20xix| - 20x 2 X2 + Uxix 2 , - 3x 2 

+10xix 2 + 5xf) + 2xi (—4x2 + x l x 2 ~ 2x 2 — xi) G(0, — x 2 , xi) 
+ (— 4xix| + 6x 2 x 2 + Ylx\x\ + x 2 — 2xix 2 — x 2 ) G(— x 2 , 0, xi) 

- (x 2 — 1) (— 6x 2 x 2 — 2x 2 + 9x2X1 + 2x2X1 — x|) H(0, xi) 

+ (— 4xix| + 6x 2 x^ — \2x\x\ + x 2 , — 1x\x 2 — x 2 ) G(— x 2 , xi)H(0, x 2 ) 
+ xi(x 2 + 1) (-4x 2 + xix 2 - 2x 2 - x i) x 

(-H(-l,x 2 )H(0,xi) -H(-l,0,x 2 ) -H(0,-l,xi)) + xi(x 2 - 1) x 
{-Ax\ + xix 2 + 2x 2 + xi) (-2H(0, Xi)H(l, x x ) - 3H(0, 0, xi) 
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+ H(0, 0, x 2 ) + 2H(0, 1, xi))} + 0(e 2 ) 



(6.34) 



, T , ,- 2£ (l-2e) JV r (l-xi)- 2g (l-x 2 )- 2e 
xk rPl =h = (s 12 ) 4 (xi+x 2 ) 2 >< 

6(xi - 1) (x 2 - 1) x 2 - 2(xi - 1) (x 2 + 2) H(0, x 2 ) x 2 
+ H(0, xi) (2x 2 H(0, x 2 ) - 2(x 2 - 1) (x 2 xi - x\ + x 2 ) ) 
+ £ jl2x 2 log(2) (xi + x 2 ) + ^7T 2 (8x1X2 - 2x2 + 4xix 2 - llx 2 + 5xi) 

+ 2x 2 G(0,-x 2 ,0,xi) + H(-l,x 2 ) (2(x 2 + 1) (x 2 xi +xi - x 2 ) H(0,xi) 

- 6(xi + 1) x 2 (x 2 + 1) ) + G(0, -x 2 , xi) (2x 2 H(0, x 2 ) 
-4 (xix| - x 2 + xi) ) + G(-X2, 0, xi) (4x 2 H(0, x 2 ) 
—2 (2x1X2 + 2x| ~~ 4xiX2 — 3X2 + Xi) ) 

+ H(-l, xi) (2(xi + 1) (x 2 - 2) x 2 H(0, x 2 ) - 6(x x + 1) x 2 (x 2 + 1) ) 
+ G(-x 2 ,xi) (12x 2 (xix 2 + 1) 
-2 (2x1X2 - 2x2 + 4xix 2 - 3x 2 + xi) H(0, x 2 ) ) 
+ 2(xi + 1) (x 2 - 2) x 2 H(-l, 0, xi) + 2(x 2 + 1) (x 2 x x + x x - x 2 ) x 
H(-l, 0, x 2 ) + 2 (x 2 + 1) (x 2 xi + xi - x 2 ) H(0, -1, xi) 
+ 2(xi + 1) (x 2 - 2) x 2 H(0, -1, x 2 ) + 6(x 2 - 1) (x 2 x x - xi + x 2 ) x 
H(0, 0, xi) - 2(x 2 - 1) (x 2 xi - xi + x 2 ) H(0, 0, x 2 ) 

- 4(x 2 - 1) (x 2 xi - xi + x 2 ) H(0, 1, xi) - 4(xi - 1) x 2 (x 2 + 2) x 
H(0,l,x 2 ) -2x 2 H(0, -1,0, xi) - 2x 2 H(0,-l,0,x 2 ) 
+ 2x 2 H(0,0,0,x 2 ) + H(0,xi) (-2H(0,-l,x 2 )x 2 + 4H(0,l,x 2 )x 2 

9 

-4(2xi - x 2 ) log(2) x 2 4(x 2 - 1) (2x 2 xi -x x - x 2 ) 

+4(x 2 - 1) (x 2 xi - x x + x 2 ) H(l, xi) + H(0, x 2 ) x 
(-2x 2 (x 2 xi - 2xi + x 2 + 2) - 4x 2 H(l, xi) - 4x 2 H(l, x 2 ) ) ) 
+ H(0, x 2 ) (4(xix 2 + 2x 2 + 3) x 2 + 4(xi - 1) (x 2 + 2) H(l, x 2 ) x 2 
-2H(0,-l,xi)x 2 - 6H(0,0,xi)x 2 

+4H (0, l,xi)x 2 + 4(2xi -x 2 ) log(2)x 2 - 2 

+ x 2 (26x 2 xi - 26xi - 26x 2 - 3C(3) + 26) } + 0(e 2 ) . (6.35) 



v( , s T , ,_ 2£ (l-2 £ ) 7V r (l-xi)- 2£ (l-x 2 )- 2g ^ 
x(-k rP3 ) =/ 5 = (-!») 2(xi+x 2 ) X 
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xi(x 2 -l) H(0,xi) -(xi-1) x 2 H(0,x 2 ) 

+ e j-^7r 2 (2x 2 xi + 5xi - 3x 2 ) - 2xi G(0, -x 2 ,xi) 

+ (2x 2 xi - xi + x 2 ) G(-x 2 ,0, xi) + (— 2x 2 xi — x\ -+ 

G(-x 2 ,xi) H(0,x 2 ) - (xi + 1) x 2 H(-l,0,xi) 

+ xi(x 2 + l) H(-l,0,x 2 ) +xi(x 2 + l) H(0,-l,xi) 

- (xi + 1) x 2 H(0,-l,x 2 ) -3xi(x 2 -l) H(0,0,xi) 
+ xi(x 2 -l) H(0,0,x 2 ) +2xi(x 2 -l) H(0,l,xi) 

- 2(xi - 1) x 2 H(0, 1, x 2 ) + H(0, xi) (2xi (x 2 - 1) 
-2xiH(l,xi) (x 2 -l) +xi(x 2 + l) H(-l,x 2 ) 

+ (xi - 1) x 2 H(0,x 2 ) - 2xix 2 log(2) ) 
+ H(0, x 2 ) (-2(xi - 1) x 2 - (xi + 1) H(-l, xi) x 2 
+ 2(xi - 1) H(l, x 2 ) x 2 + 2xix 2 log(2)) } + 0{e 2 



x 2 ) x 



(6.36) 



'6 = [8X2) 



_ 1 _ 2£ (l-2 6)jV r (l-xi)- 2£ (l-x 2 ) 
Xl + x 2 



-2s 



H(0, xi) H(0, x 2 ) + e { G(0, -x 2 , 0, xi) + G(0, -x 2 , x x ) H(0, x 2 ) 
+ 2G(-x 2 ,0,xi)H(0,x 2 ) + H(0,x 2 )(-H(0,-l,xi) -3H(0,0,xi) 



+2H(0,l,xi) 



~12 



+ H(0, xi) ( H(0, x 2 )(-2 H(l, xi) - 2 H(l, x 2 )) 



-H(0,-l,x 2 ) + 2H(0,l,x 2 ) 



7T 



+ H (0,0,0,x 2 ) 



3C(3) 



}+0(e 2 ) 



H(0,-l,0,xi)- H(0,-l,0,x 2 ) 

(6.37) 



X2 

G(-x 2 ,xi) - H(-l,xi) + H(0,x 2 ) +e{-3G(-l,-x 2 ,xi) 

- G(0, -x 2 ,xi) - G(l,-x 2 ,xi) - G(-x 2 ,-l,xi) - G(-x 2 ,0,xi) 
+ 2G(-x 2 ,-x 2 ,xi) + G(-x 2 ,xi) (-H(-l,x 2 ) + H(0,x 2 ) +2) 

- H(-l,x 2 ) H(l,xi) + H(0,x 2 )(-H(0,xi) + H(l,xi) -2H(l,x 2 ) +2) 
+ 2H(-l,-l,xi) + H(-l,0,xi) -2H(-l,l,xi) 

+ H(0,-l,xi) - H(0,-l,x 2 ) +2H(0,l,x 2 ) -3H(l,-l,xi) 
+ H(-l,xi)(-2H(0,x 2 ) +2H(l,xi) + log(2) -2) 
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+ H(l,x 1 )log(2) -^} + 0(e 2 



(6.38) 



+I{g(-x 2 ,xi) - H(-Mi)} 

- 3G(— 1, — X2,x\) + 2G(0, — x 2 ,xx) -G(l,-x 2 ,xi) - G(-x 2 , — 1, x\) 

— G(— a?2,0,zi) + 2G(-x 2 , -x 2 ,xi) 

+ G(-x 2 ,xi)(H(0,x 2 ) - H(-l,x 2 )) - H(-l,x 2 ) H(l,xi) 

+ H(0,x 2 ) H(l,xi) +2H(-l,-l,xi) + H(-l,0,xi) - 2H(-1, 1, xi) 

-2H(0,-l,xi) -3H(l,-l,xi) + H(-l,x x )(2H(l,xi) + log(2) ) 



+ H(l,xx) log(2) +0{e) 



(6.39) 



, ( ,- 1 - 2£ (l-2 £ )jV r (l-x 1 )- 2£ (l-x 2 )- 2£ 

+ I|g(-x 2 ,X!) - H(-l,x 2 ) - H(0,xi) + H(0,x 2 )} 

-4G(0,-x 2 ,xi) +2G(l,-x 2 ,xi) - G(-x 2 ,-l,xi) - G(-x 2 ,0,xi) 
+ 2G(-x 2 ,-x 2 ,xi) + G(-x 2 ,xi)(H(0,x 2 ) - H(-l,x 2 )) 
+ H(0,xi)(2H(l,xi) -3H(0,x 2 )) + H(0, x 2 ) (-2 H(l, x x ) -2H(l,x 2 )) 
+ 2H(-l,-l,x 2 ) -2H(-l,0,x 2 ) -2H(-l,l,x 2 ) + H(0,-l,xi) 
-3H(0,-l,x 2 ) +3H(0,0,xi) +3H(0,0,x 2 ) -2H(0,l,xi) +2H(0,l,x 2 ) 
-3H(l,-l,x 2 ) +3H(l,0,x 2 ) + H(-l,x 2 )(3H(0,xi) +2H(l,xi) 



7T 



+ 2 H(l, x 2 ) + log(2) ) + H(l, x 2 ) log(2) + — + 0{ e) 



(6.40) 



T { ^ 1 - 2£ (l-2e)jV r (l-x 1 )- 2£ (l-x 2 )- 2£ 
ho = { s 12) x 

Xl 



G(-x 2 ,xi) - H(-l,x 2 ) + H(0,x 2 ) + e {- G(0, -x 2 , xi) 

+ 2G(1, -x 2 ,xi) - G(-x 2 ,-l,xi) - G(-x 2 ,0,xi) + 2G(— x 2) -x 2 ,x 1 ) 
- H(0,xi) H(0,x 2 ) + G(-x 2 ,xi)(-H(-l,x 2 ) + H(0,x 2 ) +2) 
+ H(0,x 2 )(-2H(l,xi) -2H(l,x 2 ) +2) + 2 H(-l, -1, x 2 ) 
-2H(-l,0,x 2 ) -2H(-l,l,x 2 ) +2H(0,l,x 2 ) -3H(l,-l,x 2 ) 
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+ 3H(l,0,x 2 ) + H(-l,x 2 )(H(0,xi) +2H(l,xi) +2H(l,i 2 ) 
+ log(2) - 2) + H(l, x 2 ) log(2) - } + 0( e 2 



(6.41) 



'15 



_ 2 _ 2e 2- 4e (1 - 2e) N T (1 - x 1 )~ 2s (1 - x 2 ) 



-2e 



+ 2^2 + \ {h(-1,xi)+ H(-1,x 2 )- H(0,xi)- H(0,x 2 )} 
-2H(l,xi) 2 + 41og(2)H(l,x 1 ) -2H(1,x 2 ) 2 + 3G(-1,-x 2 ,xi) 

- 3G (1, -x 2 , xi) + H (0, xi) (2 H (0, x 2 ) + 2 H (1, xi)) 

+ H(0,x 2 )(3H(l,xi) + 2H(l,x 2 ))-2H(-l,-l,xi)-2H(-l,-l,x 2 ) 

- H(-l,0,xi) + 2H(-l,0,x 2 ) + 2H(-l,l,xi) + 2H(-l,l,x 2 ) 

- 2 H (0, -1, xi) - 2 H (0, -1, x 2 ) + 2 H (0, 0, xi) + 2 H (0, 0, x 2 ) 
-2H(0,l,xi) -2H(0,l,x 2 ) -2H(l,-l,xi) -2H(l,-l,x 2 ) 

- H(l,0,xi) + 2H(l,0,x 2 ) + 4H(l,l,xi) + 4H(l,l,x 2 ) 

+ H(-l,xi)(-H(-l,x 2 )+ H(0,x 2 )-2H(l,xi) + 41og(2)) 
+ H(-l,x 2 )(-2H(0,xi)-3H(l,xi)-2H(l,x 2 )+41og(2)) 

+ 4 H (1, x 2 ) log (2) + — + 0{e) J (6.42) 



1 x 



bsX\X2 



= hi = (sn)- 2 - 2 ' (- 

log (i±i) ( log ( -»-*"+' ) + log 

\l-s/\ \ 6-x i2 +2 / \6-x i2 +2 



x 



" l0g l-6,-x 12+ 2j l0g {b-x 12+ 2) + l0g(xi) 
flog f 6 r Xl2+2 ) + log f »« + *"+ 2 )) + log f \* + ^+ 2 ) x 

(_ l0g( , 2) _ 2l0g (_i^))_ 2l0g (z^^) l0g(1 _ s) 

\6 + xi 2 +2 J \ 6 + x i2 +2 J 

, / 6 s + xi 2 + 2 \ , . , / 6-6 s \ 

+ 21og h -L s + 1 -2Li 2 — 

\-6 + xi 2 +2y \6-xi 2 +2y 

/ -6-6s \ / 6s-6 \ / 6-6s \ 

+ 12 V-6 + xi 2 + 2y ~ 12 V-6 + xi 2 + 2/ + 12 V& + X12 + 2 / 
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where 



2 Li, 



2 



+ 2Li 2 



6 + 6s 
6 + X12 +2 
Xl2 — 6s 

6 + X12 
6s + Xi 2 
6 + X12 



+ 2Li 2 



-bs + xu - 2 



2 Li. 



6 + x 12 - 2 
6 s + X12 - 2 



2 Lis 



X12 - 6s" 



+ 0( 



6 + X12 - 2 



+ 2Li 2 



»12 
6s + X12 
X\2 - 6 



(6.43) 



X\2 = Xi + X 2 

6 = 



12 



+ 4. 



Q 3 = xi (xf 2 + 4) - 4xi 2 



The master In was calculated using differential equations and is written in terms of log- 
arithms and dilogorithms of complicated arguments instead of GHPL's and HPL's. The 
coefficients of the differential equations for In involve a new denominator, cubic in x\ and 
quadratic in X2, 

Q 3 = xi (xi + x 2 ) 2 - 4 x 2 (6.44) 

which does not occur in the differential equations for any of the other masters presented 
in this paper. Because of this denominator, In has this more complicated form. 

The differential equations for In are solved by first considering the corresponding 
homogeneous equations. A solution of these equations, for e = 0, is given by 

I\r = ^Jfy (6-45) 



Note that it is singular along the curve Q3 = 0, part of which is inside the physical region 
< x\ < 1, < X2 < 1. On the other hand, the full solution for In, which we write as 

In = Cnl\2 m (6-46) 

cannot be singular inside the physical region. Therefore, the coefficient Cn must vanish 
along the curve Q3 = 0. In order to solve the differential equations for In-, it is convenient 
to consider them keeping X12 = x\ + X2 fixed, so that Q3 = xi(x\ 2 + 4) — 4xi2 becomes 
linear in x\. Using the boundary condition that Cn must vanish at the point x\ = x\ = 
4xi2/ (xf 2 + 4), where Q3 = 0, the solution of the differential equation of Cn in x\ can 
then be written as a one-dimensional integral: 

Cn - 

Jx\ V^iV^H^ +4) - 4xi2 

log(x'i) 
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+ (log(l + X[) - log(x 12 )) + log(xi 2 - X' x ) 



.47) 



The integration can be performed analytically in terms of logarithms and dilogarithms. 
The final result for /14 is shown in eq. (6.43). 

We have presented the result for In in a form where all (di) logarithms are real in the 
region where Q3 > and s is real. By construction, J14 has no singularity at Q3 = 0. 
For Q3 < 0, s is imaginary and the arguments of the (di) logarithms that depend on s 
become complex. Nevertheless, I14 remains real in this region because the factor C14 is 
now imaginary. 

Finally, it is worth noting that I14 is only needed up to finite order in the hard region. It 
is required up to order e in the collinear x\ region where it has been calculated separately 
using differential equations and direct integration. 

b) The collinear regions 

The real-real master integrals were calculated in the collinear regions by deriving a hyper- 
geometric integral representation starting from their definition as phase space integrals. 
This allows us to give the results in closed form. We also have derived these results as 
an expansion in e using the differential equations method, which provides an important 
check on our results. We list below the results of the masters contributing to the antennae 
£>4, £4, %\ in these regions. 

• collinear x\ 

In this region, only h, I 2 , I7 and Is require to be expanded up to C(e 3 ) since they contribute 
to E12 which starts its power expansion at C(e -3 ). 



k k -pi =h = (s 12 y- 2£ N r 2-^(l-x i y-^(l-x 2 y- 2£ x 2 s x (6.48) 



2 



1 + 2j r(2-2e) 

I 

^Ql^ =h = {si2)- 2e N T 2~ £ (1 - Xl ) l - 2s (1 - x 2 y- 2£ x 2 s (1 + x 2 y 1+e x (6.49) 
i 

r(i-g) 2 

T(2-2e) 




>^x(-(Vpi)(VP3)) =h = (s 12 ) 1 - 2£ N r 2-^(l-x 1 ) 1 -^ x 



(6.50) 



(l-x 2 ) 1 - 2£ x^ e (l+x 2 )- 1+e X 
(1 - 3x 2 + 2xj 2F1 (1,1 - e;2 - 2e; 1 - x 2 )) 



r(l-e) 




T(2-2e) 

■/j^xkj-pi =h = -(s 12 y 2£ N r 2- 1 - £ (l-x 1 ) 1 - 2£ (l-x 2 ) 1 - 2£ X2 e x (6.51) 

T(l-e) 2 



(1 + x 2 )~ 1+£ (-1 + x 2 2 F X (1, 1 - e; 2 - 2e; 1 - x 2 )) 



r(2-2e) 




X x (-fcj • p 3 ) =h = -{s 12 y 2£ N r 2' l - £ x 



(6.52) 




(i- Xl y- 2£ (i-x 2 ) L - 2£ x 2 £ (i + x 2 y 

(-l + x 2 2 F 1 (l,l-e;2-2e;l-x 2 )) 



h = (s 12 )- L - 2£ N T 2- £ (1 - x x y- 2£ (1 - x,) 1 "^ x 2 £ x 



-l+e 



r(i-e) 2 

r(2-2e) 



(6.53) 



(l + x 2 ) 



-l+e 



r(i-e) 2 



2 Fi(l,l-e;2-2e;l-x 2 ) 





r(2-2e) 

/7 = -(si 2 )- 1 - 2e iVr 2-^ £ (1 - xx)~ 2e (1 - x 2 ) l ~ 2£ x 2 x - £ x 

(1+ , 2) . 2Fl ( M _ E;2 _ 2s; l-) ^ 

/ 8 = M' 1 ^ N r 2- l - £ (1 - Xl y 2£ (1 - x,) 1 - 25 x 2 l - £ x 



/ -l + x 2 \ 1 r(l-e) 2 



(6.54) 



(6.55) 



'14 



(s l2 )- 2 - 2£ N r 2^- £ (1 - Xl )- 2£ (1 - x 2 ) x ^ £ x 2 L - £ x 



2e 



vl-2e -1-e 



< 



{l + x 2 f 2 F x ( l,l-e;2-2e;^-^ 2 - ) 2 F X (1, 1 - e; 2 - 2e; 1 - x 2 ) 



(6.56) 

r(i-g) 2 

r(2-2e) 



collinear x 2 



In this region, only the following four masters (all of them are needed up to 0(e 2 )) con- 
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tribute: 



y^X^)\ xkk ' Pl =h = (su) 1 - 2£ N r 2~ 3 - £ (l-x 1 ) 1 - 2s x 



(1 - x 2 f- 2e x7 £ (1 + Xl )- 2+£ (1 + 3x0 



h = {si2)- 2e N V 2- £ (i - Xl ) l - 2e (i - X2 y- 2£ xr e x 



\l-2e „-e 



Jl + xi)" 



r(i - e) 2 

r(2-2e) 




(6.57) 



r(i-e) 2 

r(2-2e) 



/ 7 = -( Sl2 )- 1 - 2e Ar r 2^(l-x 1 )- 7e (l-x 2 ) i -^xr £ x (6.59) 



2e 



\l-2e -e 



(6.58) 



(1 + xi)" 



-i +£ r(i - e) 2 



r(2-2e) 




'10 



(s 12 )- 1 - 2£ n t 2-^ e (i - Xl y-' ze (i - x 2 )-^ e xr 1 ^ x 



1~2 £ 



-2e „-l-e 



(1 + Xi) f 



r(i-g) 2 

r(2-2e) 



Fx l,l-e;2-2e; 



l-a?l 



(6.60) 



c) The soft region 



As discussed previously, the master integrals in the soft region are expected to be needed at 
most two orders in e higher than in the hard region. It turns out here that the masters are 
only needed up to C(e 3 ) at most in the soft region. The reason behind this is the absence 
of double soft gluon configurations in the antennae B, E and H. Furthermore, as can be 
seen from Table 1, only B\ 2 will have a contribution from the soft region at all. This is 
due to the fact that only B\ 2 allows for a singular double soft configuration, a soft q — q 
pair. The master integrals which are required in this case up to 0(e 3 ) are either I 2 or 
depending on the basis choice for the masters. The others are only needed up to 0{e 2 ). 

In the soft region, the masters are calculated by a direct evaluation of the phase space 
integrals. In this case, the results are products of gamma functions and therefore can be 
presented in closed form. 
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x k k - Pl =h = {s l2 y- ie N T (i - Xl y- 2£ (i - X2 ) 



r(i 



\2-2e 




8T 2 - 2e 



/ 2 = (si 2 )- 2e iV r (l-x 1 ) 1 - :4E (l-x 2 



1-2 E , 



a-2e 



r(i 



2r(2 -2s) 



h = -(s 12 )- 1 - 2e N r (1 - Xl )- 2£ (1 - X2 y 



(6.61) 



(6.62) 



(6.63) 




r(i-g) 2 

2T(2 - 2s) 
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(»i2)" 1 ~ 2e JVr (1 - si) 1 " 26 (1 - x 2 y Z£ x 



1-2 E 



-2e 



r(i-^ 



(6.64) 



2(1 - 2e)r(l - 2e) 

6.2.2 Integrated antennae 

In this section, we present the results for the integrated crossings of three final-final four 
parton antenna functions B®(q,q' ,q' ,q), E®(q,q' , <f , g) and H®(q,q,q' ,q') defined in [19]. 
Because the results are too long to present fully here, we show explicitly only the pole 
terms (all pole terms higher than l/e 2 , and only the highest pole term if it is 1/e). The 
dimensionful prefactor (si 2 )~ 2e of every integrated antenna is omitted below. The complete 
results are included in a Mathematica file appended to the source file of the manuscript. 

The integrated form of the initial-initial antenna functions Hi 2 and H\% are ob- 
tained by first crossing a pair of identical or non-identical quarks in the final-final antenna 
H® (lq, 2q, 3q' , Aq') to the initial state, and then calculating the phase space integrals as 
described in the previous sections. They take the following form: 



n 13 



'\X\X 2 + 1) {{x\ + 



cj - 4) x\ + l) x\ + xQ 



3 (xi + x 2 ) 4 
2xi + 2) (xl - 2x 2 + 2) 
Ax\x 2 



} + o(i) 



(6.65) 
(6.66) 



For the E^flq, 2q' , 3q' , 4g) antenna, there are four independent expressions, obtained by 
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crossing qq' , qg, q'q' or q'g to the initial state, we list the pole terms of the integrated ones: 
1 r(x 2 -2x 2 + 2) 6{1 -xi) 



£ 



12 



I 4x 2 J 



+ F { " 2(X 1 + 1)X 2 (X 1 +X 2 )2 (" 2X ' X 2 " 2*1X2* " X^ + 2X!X^ 

+2x2 + x i x i ~ 2xix 2 - 2x 2 - 2x 2 x 2 + 2x\ - 4xix 2 + 4xi) 

(x| - 2x 2 + 2) Pq(xi) _ / (x 2 -2x 2 + 2) ff(-l,x 2 ) 

2x 2 1 Xlj ^ 2x 2 

(3x2 _ q X2 + g ) ^ ^2 _ 2x2 + 2 ) 



8x 2 2x 2 
81og(2)xl-7xi-161og(2)x 2 + 24x 2 + 161og(2) -20^ /1 



16x 2 



! )} 



1_ _ xi) ( (x 2 + l)H(0,x 2 ) _ (x 2 -l) (4x| + 7x 2 + 4) \ i + ^ ^1 



£u = ^ 5(1 - xx) v V - — ^rrf f +° - K 6 - 68 ) 



1 rx 1 x 2 (x 1 x 2 + l) 2 (x 2 + x\ - 2) 



^24 = - (- 7 », 1 \(x 2 ~ 1) (I2x? + 36x 2 x^ - 6xf + 41a%x[ - 19x 2 x^ (6.70) 

£ L 6xJ(xi + X2) I 

-x\ + 19x^xf - 29x^xf - 5x 2 xf + 4x^xf - 20x^xf - 26x|xf + 15x 2 xf 
-Qx\x\ - SOx^xf + 13x\x\ + 10x 2 x? - 2x\ - \2x\x\ + 86x^xf + 50x^x? 
-34x 2 x? + 32x\x\ + 62x 2 i x 2 - \\%x\x\ + 24x^x x - V&2x\x x - 48x^) 



(2xf + x 2 xj - x\ + 2x 2 xf - x 2 xf - 8X2X 2 + 2x 2 xi + 8x3) 

x\ 

(G(-x 2 , xi) - H(-l, xi)) + x\ (2x\ + X!X 2 - 1) H(0, x 2 )} + O (1) . 

Finally, for the B®(lq, 3q', 4q', 2q) antenna, we have three independent crossings, obtained 
by crossing either the primary quarks qq, the secondary ones q'q', or a combination of 
primary and secondary quarks to the initial state. Below we show the pole terms of B\ 2 
and i?i3, the antenna -B34 is completely finite. 

6i2 = _ ji | *(i-*,H(i-x 2 ) | (6 . 71) 



4{*d-x 1 )(-i4p + ^„(x 2 )-l (S (i-^ 

+ (i B „(x 1 ,-l^l)*(l-x 2 )} +0 (i), 
Ba _ -Li(l - { (1 - I2>(4 2 l +7X2+4) + ^(0,x 2 )} (6.72) 



,2 
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7. Conclusions 



Within the antenna subtraction formalism, allowing the calculation of higher order QCD 
corrections to jet observables, subtraction terms are constructed from antenna functions. 
Those functions describe all unresolved radiation between a pair of hard radiator partons. 
At NNLO, this formalism has been fully developed and applied so far only for colour- 
less initial states. In this paper, we have focussed on the extension of this formalism to 
evaluate NNLO corrections to jet observables at hadron colliders and concentrated on the 
construction of subtraction terms for the double real radiation contributions. More pre- 
cisely, we have considered the subtraction terms needed to account for the radiation of two 
colour-connected unresolved partons off two initial state partons. For these subtraction 
terms, four-parton tree-level initial-initial antenna functions are required in unintegrated 
and integrated form. The integration over the phase space associated with two unresolved 
partons has to be performed analytically. 

In this paper, we have given a catalogue of all non-identical four-parton initial-initial 
antenna functions. Furthermore, after applying standard reduction techniques, we found 
that 32 master integrals are necessary to obtain their integrated form. As a step towards 
the integration of the full set of integrated initial-initial antenna functions, in this paper 
we have focussed on the initial-initial antennae obtained from crossing two partons in the 
final-final antenna functions characterised by the presence of two quark flavours. After 
reduction, 12 masters were required to obtain those. We presented the decomposition of 
the calculation according to four phase space regions: hard, collinear and soft and we gave 
the master integrals in these regions. Finally, we presented the results for those integrated 
initial-initial four parton antenna functions themselves. Since the results are lengthy, we 
have shown only the leading pole terms in the manuscript and have attached the complete 
results as a Mathematica file. 
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